Low-energy instability of flexural phonons in graphene 
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We study the effect exerted by the electrons on the flexural phonons in graphene, accounting for 
the attractive interaction created by the exchange of electron-hole excitations. Combining the self- 
consistent computation of the phonon self-energy with renormalization group methods, we show that 
graphene has two different phases corresponding to soft and strong renormalization of the bending 
rigidity in the long-wavelength limit. In the first case, the system may have an intermediate scale 
in which the phonon dispersion is softened, but it manages finally to become increasingly rigid over 
large distance scales. The strongly renormalized phase is closer however to critical behavior, with 
an effective rigidity that becomes pinned in practice at very small values, implying a very large 
susceptibility for the development of a condensate of the flexural phonon fleld. 



In recent years, the feasibility of fabricating graphene 
layers a single atom thick[T]-(3] has led to a great deal of 
activity, because of their novel electronic properties and 
potential applications [4]. Beyond its remarkable elec- 
tronic transport attributes, it is becoming increasingly 
clear that graphene's mechanical and electromechanical 
properties are just as remarkable. Despite it's extreme 
thinness, it behaves as a stable metallic membrane with 
finite elastic moduli and rigidity. From this point of view, 
one of the most intriguing features of graphene is its 
tendency to develop ripples with a characteristic length 
scale. In exfoliated graphene, these ripples are correlated 
to some extent with the irregularities of the substrate 
on which the graphene sheet is deposited. But there is 
also evidence that they arise in part as an effect intrinsic 
to the two-dimensional membrane. When graphene was 
first fabricated, it came as a surprise that a purely two- 
dimensional material could exist, without being destabi- 
lized by large thermal or quantum fluctuations. From 
a theoretical point of view, it has been shown however 
that the coupling of flexural phonons to the in-plane vi- 
brations in a crystalline membrane provides a mechanism 
to stabilize the two-dimensional system, leading to a fi- 
nite critical temperature for the crumpling transition [7]. 
The ripples are usually viewed then as a sort of imprint 
of the large fluctuations that should be present in the 
low-dimensional material [8]. 

The ripples are expected to have a significant impact 
on the electronic transport in graphene [3]. More gener- 
ally, it has been shown that the scattering with flexural 
phonons may lead to a reduction of the lifetime of elec- 
tron quasiparticles in graDhene|10). On the other hand, 
the effect exerted by the electrons on the out-of-plane 
vibrations of a metallic membrane has been much less 
studied. There have been partial approaches in which 
the coupling of the charge density to the flexural phonons 
has been considered |l H - [T3] . but a zero temperature scal- 
ing analysis characterizing the possible critical behavior 
in the system has not been undertaken until now. 

In this paper we investigate the many-body effects on 
flexural phonons in graphene, in a theory where these are 
coupled to the electron charge density as well as to the 



in-plane vibrations. These are known to give rise to a 
repulsive self-interaction between flexural phonons. We 
show that the exchange of electron-hole excitations leads 
to an interaction with the opposite character, although it 
plays in general the role of an irrelevant perturbation in 
the system. We will analyze the conditions in which this 
attraction may have a significant effect, finding that it 
leads in that case to a softening of the acoustic branch of 
out-of-plane phonons in graphene. The effective bending 
rigidity can be then strongly suppressed beyond a cer- 
tain length scale, making the system very susceptible to 
the development of a condensate (nonvanishing average 
value) for the out-of-plane phonon field. 

We characterize the shape of the graphene sheet by 
the vector field u = {ui,U2,h), where ui,U2 represent 
the in-plane displacement with respect to the equilibrium 
position and h is the out-of-plane shift. The deformation 
of the membrane is analyzed in terms of the strain tensor 
Uij, given by 

= ^{diUj + djUi + dihdjh) (1) 

Thus, the action of the system has a term governing the 
dispersion of flexural phonons and depending on the mass 
density p and the rigidity k 

5ft.ce = \jdt (fx{p{dthf - K{V^hf) (2) 

while the dynamics of the in-plane phonons is controlled 
by a term of the form 

Sn = ljdt (fx{p{dtu,f ~ 2pTr u% - \{Tr u,jf) (3) 

where p and A are the Lame coefficients. 

The term 5*1, contains a quartic self-interaction for the 
out-of-plane phonon field h. It is well-known that this in- 
teraction has to be taken into account for the description 
of stable crystalline membranes, as the soft dispersion in 
([2]) leads to infrared divergences that cannot be tamed in 
the free theory]^. The novel effect arising in graphene 
is that the electronic charge couples to the displacement 



of the sheet, in such a way that it gives rise to another 
source of self-interaction between the phonon fields. We 
will consider the effect of electrons from the tt bands of 
graphene, represented by the field ^(x). Accounting for 
the fact that the coupling must take place through the 



trace of the strain tensor it. 
term 



we add to the action a new 



5*0 



-ph 



(4) 



The electron-phonon coupling g has dimensions of energy 
and is directly related to the deformation potential of the 
graphene lattice, which has a microscopic origin. 

Assuming that the electron correlations are weak in 
graphene, we may represent the dynamics of the 'i!{x) 
field by a quadratic hamiltonian, matching at low ener- 
gies the usual expression for the Dirac fermion fields. The 
electron degrees of freedom can be integrated out at this 
point, leading from (4]) to a new phonon contribution to 
the action, to be added to ^ and ([s]), 

Sn' = -^g'^ J duqd^qx(Si,^^q)(Tru,j){Tru^.j) (5) 

where Wg) is the density-density correlator depend- 
ing on frequency and momentum variables. Terms that 
are higher order in Tr Uij can be generated after inte- 
gration of the electron fields, but they turn out to be 
more irrelevant perturbations in the approach developed 
in what follows. 

Since the energy scale of electronic excitations (for a 
given momentum) are much larger than the phonon en- 
ergies, we can safely take the limit a; = in the charge 
susceptibility. In the Dirac theory, we have in particular 
x(q, 0) = — |q|/vi?. We observe that the new contribu- 
tion (pi tends to cancel the term with the LamE coef- 
ficient A in ([3|. After following the usual procedure to 
integrate the in-plane phonons Ui [7], we get the total 
action for the fiexural phonons 



dujq d^q K{q)uij{q,u>g)uij{-q, -ujq) (6) 



where Uij{q,Uq) is the Fourier transform of 
{l/2)Pijdihdjh, Pij being the transverse projector, 
and the quartic coupling is (in the Dirac theory) 



K{q)^ 



Afi{fi + X-g'^\q\/vF) 
2fi + X~g-^\q\/vF 



(7) 



Negative values of the coupling function K{q) may lead 
to an instability of the graphene sheet, as they represent 
an attractive interaction between the flexural phonons. 
The phonon self-energy and the interaction vertices are 
actually affected by logarithmic divergences in the energy 
cutoff, which makes it necessary to keep track of the low- 
energy dependence of the couplings. In order to study 
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FIG. 1. First order corrections to the self-energy of flexural 
phonons (represented by a wiggly line) arising from (a) the 
four-phonon interaction Kq and (b) the exchange of electron- 
hole excitations (represented by the bubble with arrow lines) . 

the competition between positive and negative couplings 
in ([?]), we may approximate this coupling function by the 
constant term Kg = 4/i(/i + A) / (2/i+ A) and the dominant 
powers of |q| 



Kiq) ^Ko-Gi \q\ - q' 



(8) 



Neglecting higher-order terms in |q| amounts to disre- 
garding the exchange of an increasing number of electron- 
hole pairs in the in-plane phonon propagator, which is 
justified as long as the pole arising from the denomina- 
tor of ([t]) is not included within the first Brillouin zone 
of a real graphene system. This is consistent with an 
estimate for the on-site deformation potential for the in- 
plane phonons g^^ « 8 eV [14j, that is the coupling ap- 
pearing in the denominator of Eq. ([t]) . 

The different interactions tend to modify the low- 
energy behavior of the rigidity k in opposite directions 
depending on their attractive or repulsive character. This 
can be seen through the corrections to the phonon self- 
energy S(p,Wp) represented in Fig. [T] Computing the 
full propagator of the flexural phonons from the expres- 
sion D~^{Tp,u}p) = ujp — {ko/p)p'^ — I](p, Wp), we find the 
renormalized bending rigidity k(p) given by the equation 



«(P) 



P 
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d(t) 



(9) 



sin4((/)) 



/\o-Gi|p-q|-G'2|p-q|' 



p2|p. 



A self-consistent treatment of the many-body theory re- 



quires that the rigidity arising in the integrand of ( 10 ) 



from the phonon propagator must be also taken as k(pJ 
This leads to an integral equation, that can be solved to 
obtain the effective momentum dependence of the renor- 
malized parameter. 

We have represented in Fig. [2ja) the result of solving 
Eq. ( 10 ) in different regimes of the coupling g. We ob- 
serve that, for a sufficiently large value gc, the function 
k(p) vanishes at a certain momentum, bouncing back for 
smaller values of |p|. In this picture, gc plays the role of 
critical coupling, as larger values of g would actually lead 
to negative values of the bending rigidity, implying the 
absence of a self-consistent solution. In that case, how- 
ever, a more refined analysis has to be applied in order 
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FIG. 2. (a) Effective momentum dependence of the bending 
rigidity obtained from the solution of Eq. ( 10 1 for values of Gi 
and 6*2 corresponding (from top to bottom) to a deformation 
potential g of flexural phonons equal to 0, 24 and 26.3 eV. 
(b) Scaling of k in the renormalization group approach, for 
bare values of the deformation potential equal (from top to 
bottom) to 0, 22 and 29 eV. In both approaches, we have 
taken the bare value of a^Ko equal to 18 eV. 



to capture the low-energy scaling of the couplings Kq, 
Gi and G2, which becomes quite significant when the 
effective bending rigidity gets very small. 

We apply then a Wilsonian renormalization group 
method starting from a high-energy cutoff Ec = yj kq / pq^ 
for some large momentum (jc, and performing a progres- 
sive integration of energy shells to approach the low- 
energy regime. The bare coupling Kq is corrected to 
lowest order by the one-loop diagram made of the ex- 
change of two flexural phonons, which shows a logarith- 
mic dependence on the cutoff q^. This can be absorbed 
into a suitable renormalization of the effective coupling. 
Under a differential change of the energy scale, we find 
the scaling with the running cutoff 5 — > 



dq 



3 {Ko/p^y 
647r (k//9)3/2 



(10) 



The couplings Gi and G2 correspond to effec- 
tive phonon interactions mediated by the exchange of 
electron-hole excitations, and they also undergo a log- 
arithmic renormalization from the exchange of flexural 
phonons. However, the main source of scaling for these 
couplings arises at the classical level, as the respective 
terms in the effective action pick up anomalous factors 
^/y/s and 1/s when shrinking the energy scale from Ec to 
Eel s. We find then the respective anomalous dimensions 
A = —1/2 and —1, which lead to the scaling equations 
for 71 = 1 , 2 



n G„ 



3 {Ko/p^)Gn 



dq ' 327r (k/p)3/2 

The main scaling dependence is given by the behavior 



(11) 



where the couplings G„ are subject to a purely logarith- 
mic scaling in the same fashion as Kq. 

A scaling equation can be also written for the rigidity 
K, as the self-energy diagrams in Fig. [T] actually display 
a logarithmic dependence on the cutoff qc- Taking into 
account the irrelevant character of the couplings G„, we 
find that under a variatioir of the momentum scale q 



Ok/ p 
' dq 



3 Kq/p" 3 



Eq. ( 13 1 has to be solved then together with ( 10 ) and 



(111 in order to obtain the behavior of the effective rigid- 
ity K. We note that, in real graphene samples, the bare 
values of the couplings Kq and G„ place in general the 
system in the weak-coupling regime. On the one hand, 
the energy scale of Kq is given by a?Ap{p + X) / {2pi-\-\) « 
20 eV jT5], where a = 1.4 A is the interatomic spacing. 
The coupling from the electronic sector is subject to the 
uncertainty in the estimates of the deformation potential 
for out-of-plane vibrations, leading to a wide range of 
values with g ~ 20 — 30 eV [TB]. On the other hand, we 
have the energy scale c? p^(k.I p)^!"^ ^ which has a nominal 
value of « 40 eV in graphene. When nj p becomes very 
small, however, the rate of decrease of Kq^ G\ and Gi is 
greatly increased, leading to significant differences with 
respect to the behavior of the effective rigidity obtained 
from the self-consistent equation (10). 



The behavior of k obtained from the scaling equations 
is shown in Fig. [2|b) for different bare values of g. For 
small values of this coupling, we see that the effective 
rigidity has a steadily increasing trend as q — > 0. There 
is however a regime of large 17, corresponding to values 
of the deformation potential g ~ 25 — 30 eV, in which 
K is significantly reduced due to the electron-phonon in- 
teraction. At some point, k may be even suppressed by 
several orders of magnitude, but without vanishing as in 
the self-consistent approach shown above. Here we find 
instead that k always bounces back from zero, though 
it may become trapped in a range of very small values 
as the effective couplings Kq^G\ and G2 are themselves 
quickly renormalized very close to zero. 

We find therefore that graphene has two different 
phases corresponding to soft and strong renormalization 
of K as (7 — )• 0. In the first case, the system may have 
an intermediate scale in which the effective k is reduced, 
but it manages finally to become increasingly rigid over 
large distance scales. In the strongly renormalized phase, 
however, the effective rigidity becomes pinned in practice 
at very small values below a certain momentum scale. 

Although the self-consistent solutions of Eq. ( 10 1 miss 



G„((z) = g"G„(9) 



(12) 



the renormalization of the couplings at very small k, we 
can still use them to draw the boundary between the 
phases with soft and strong renormalization of the bend- 
ing rigidity. That line can be formed with the values of 
critical coupling gc beyond which the integral equation 
would lead to negative values of k. This forbidden region 
corresponds in the renormalization group approach to the 
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FIG. 3. Diagram showing the phases corresponding to a rigid 
long-distance regime of graphene (lower region) and pinning of 
the renormalized bending rigidity at very small values (upper 
region), in the space of nominal values of the rigidity kq and 
deformation potential g. The value of a^Ko was taken as 18 
eV. 



regime where the bending rigidity becomes trapped at 
very small values. The resulting phase diagram in the 
space of nominal values of the rigidity kq and deforma- 
tion potential g is shown in Fig. [3j 



The phase with strong renormalization of the bending 
rigidity can be put in correspondence with the appear- 
ance of ripples in graphene, as very small values of k 
imply a very large susceptibility for the development of 
a condensate (nonvanishing average value) of the flexu- 
ral phonon field h. It is conceivable that the different 
degree of Coulomb screening may control which of the 
two phases is found in a real graphene sample, as the 
Coulomb interaction leads to a slight reduction in the 
charge susceptibility. Thus, graphene can be more prone 
to develop ripples or to buckle when surrounded by me- 
dia with high dielectric constant, for which the values of 
x(q, 0) are maximized. 
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